Magnetic properties of superconducting multifilamentary tapes in perpendicular field. 

I: Model and vertical stacks 
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Current and field profiles, and magnetization and AC losses are calculated for an array of infinitely 
long superconducting tapes arranged vertically in a perpendicularly applied magnetic field. Calcu- 
lations are based on the critical state model. The finite thickness of the tapes and the effects of 
demagnetizing fields are considered. The influence of the magnetic coupling of the filaments in the 
magnetic properties of the arrays are studied. The general model can be applied to an arbitrary 
arrangement of tapes as long as there is reflection symmetry with respect to the vertical central 
plane. 
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I. INTRODUCTION 

The magnetic behavior of hard superconductors de- 
pends not only on their intrinsic properties, such as the 
critical-current density J c , but on their geometry. For 
example, the magnetization and ac losses for a supercon- 
ducting thin film are very different depending on whether 
the magnetic field is applied perpendicular or parallel to 
the thin dimension. Recent advances in the modeling of 
magnetic properties of superconductors (SCs) have al- 
lowed to understand their behavior in a large variety of 
geometries. The most common framework used to ad- 
dress this problem has been the critical state-model [Q, 
that assumes that currents circulating in the SCs flow 
with a constant density J c , later extended to currents 
depending only upon the local magnetic field Hi Q . The 
original model was solved in the parallel geometry, that is 
for applied fields H a along an infinite dimension for slabs 
and cylinders , because in these geometries the prob- 
lem of the demagnetizing effects was not present. A fur- 
ther step was presented when the CSM was extended to 
the case of very thin strips Q-j^] and disks , for which 
important demagnetizing fields were involved. More re- 
cently, the more general case of critical state in samples 
with finite thickness, such as strips Jlcfl and cylinders 
[p~T|— [l3[| was solved by numerical models. 

However, there is a particular arrangement that has 
not been solved so far in spite of its practical impor- 
tance. It is the case of a set of superconductors with 
finite thickness arranged in an array form. The impor- 
tance of this case is not only academic but practical, 
since this geometry is often a good system for model- 
ing the filamentary structure of actual superconducting 
tapes. The recent advances in the technology of super- 
conducting Ag/Bi2223 multifilamentary tapes increase 
the importance of having a theory describing the mag- 
netic properties of such arrays. 

Although the problem of studying the current distribu- 
tion and magnetization in array of superconducting tapes 
has not been systematically solved, there have been sig- 
nificant works offering partial solutions. Fabbricatore et 
al jl4j presented a comprehensive analysis of the Meissner 



state in arrays of strip lines arranged vertically (z stack of 
strips), horizontally (x arrays) and in the form of a ma- 
trix (xz array) and compared their results with actual 
measurements on multifilamentary tapes. Their numer- 
ical procedure, however, was not adequate to study the 
more general case of bulk current penetration. Mawatari 
Jl5| studied not only the Meissner state but also the crit- 
ical state for the case of an infinite set of periodically 
arranged superconducting-strip lines, in the limit that 
the strip lines were infinitely thin. Mawatari and Clem 
Jl6| studied the penetration of magnetic flux into current- 
carrying (infinitely thin) strips lines with slits in the ab- 
sence of applied magnetic field. All the existing models 
assume either arrays of infinitely thin strips in the criti- 
cal state or arrays of strips with finite thickness but only 
in the Meissner state. The only exception we know is the 
recent work by Tebano et al [ jl7| in which preliminary re- 
sults on the current penetration and magnetization were 
calculated for some realistic arrays based on the proce- 
dure developed by Brandt JlO| | . 

A key issue in the study of superconducting tapes 
is how currents circulate within filaments. There are 
two important cases concerning this point, depending on 
whether current in each filament is restricted to go and 
return through the same filament or if instead there is no 
such restriction. In principle, both cases have practical 
and theoretical interest, as we discuss in the second paper 
of this series . The magnetic response and AC losses 
will strongly depend on this factor so it is imperative to 
study both cases separately and compare them. 

In this series of papers we study the current and field 
penetrations, magnetization and AC losses of arrays of 
superconducting strips of finite thickness. In the first pa- 
per we present the model and its application to the case 
of an array of a finite number of infinitely long strips 
of finite thickness arranged vertically (z stack of finite 
strips) with a perpendicular applied field. This geometry 
is studied first because it is independent of the connec- 
tion type, since because of the system symmetry, current 
always go and return through each strip. In the second 
paper we will study the cases of horizontal (x) arrays 
and matrices (xz array), for which different behaviors 
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arise depending on the connections. In all cases, we will 
concentrate our study in tapes composed of strips with 
high aspect ratio since this is the case most often met 
in practice, although our model is applicable to arrays of 
strips with arbitrary thickness. However, this is differ- 
ent from the approximation of considering infinitely thin 
strips as in Jlq] , since we take into account the different 
current penetration across the superconductor thickness. 

The present paper is structured as follows. In section 
II we introduce the calculation model. Current and field 
profiles are calculated and discussed in Section III. The 
results of magnetization and ac losses are discussed in 
Section V and VI, respectively. In Section VII we present 
the main conclusions of this work. In two appendices we 
present analytical expressions for both the full penetra- 
tion field of the arrays and the inductances used in the 
model calculations. 



II. MODEL 

The model we present here is suitable for any supercon- 
ductor geometry with translational symmetry along the 
y axis and mirror symmetry to the vertical plane. How- 
ever, we will focus on z-stacks, x-arrays, and xz-matrices 
made up of identical strips infinitely long in the y direc- 
tion. The horizontal and vertical strips dimensions are 
2a and 2b, respectively. We consider a uniform applied 
field H a in the z direction. Results for z-stacks are dis- 
cussed in the present paper, whereas those for x-arrays 
and ccz-matrices are presented in the second paper of this 
series @. 

Our numerical model is based on minimizing the mag- 
netic energy of the current distribution after each applied 
field variation. The model assumes there is no equilib- 
rium magnetization in the superconductor and, in the 
present series of papers, we further assume there is no 
field dependence of J c for simplicity. The energy and 
flux minimization in the critical state was previously dis- 
cussed by Badia jL9) and Chaddah et al The 
details of the general numerical model can be found in 
p2|,fL3 23 2m. The approach has been successfully ap- 



plied to describe the experimental features observed in 
the initial magnetization slope (2^] as well as the whole 
magnetization loop and levitation force [|1||3| of super- 
conducting cylinders. We now outline the main charac- 
teristics of the model. 

For cylindrical geometry j2^JI^,^], the superconduct- 
ing region was divided into a certain number of elements 
of the shape of rings with rectangular cross-section, which 
formed circular closed circuits. For a given current dis- 
tribution in the cylinder, the energy variation of setting 
a new current in a certain element was calculated. The 
minimum energy variation method consists in, given an 
applied field H a , looking for the circuit which lowers the 
most the magnetic energy and set there a new current of 
magnitude J c , and then repeat the procedure until set- 



ting any new current does not reduce the energy. The 
method as explained above is used to find the initial 
magnetization curve Mi(H a ). Provided that J c is field 
independent, the reverse curve can be found using that 



M rev (# a ) = Mi(H m ) - 2M l ({H m - H & )/2), (1) 



and the returning curve using 



M ret (H a ) = -M rev (-tf a ), 



(2) 



where H m is the maximum applied field in the loop. 

In the present case, each strip is divided into a set of 
2n x x 2n z elements with cross-section (Ax)(Az) and an 
infinite length in y direction, as shown in Fig. 1. The di- 
mensions of each element are Ax — a/n x and Ay = b/n y . 
We consider that the current density is uniform within 
the elements and flows through the whole element section 
and not only through a linear circuit, as in p2] , p^| , p3[ . 

When the distribution and orientation of the strips in 
the set are symmetrical to the zy plane, as is the case 
for a z-stack, in the presence of a uniform applied field 
the induced current front will also be symmetrical to the 
zy plane. Thus, we can consider that the pair of ele- 
ments centered at (x, z) and (— x, z) form circuits that are 
closed at infinity. This grouping in pairs, forming closed 
circuits, allows for the analytical calculation of self and 
mutual inductances per unit length of the circuits with 
finite cross-section (see the Appendix B for inductances 
derivations and formulae). The pairs, or circuits, are la- 
belled using the subscript % from 1 to N = 2n x n z n,f, 
being n f the number of strips of the set and N the total 
number of elements in the x > portion of the set of 
strips. 

Once the analytical expressions for the inductances are 
obtained, the energy of the i circuit can be calculated as 



N 

£ 



Mijljli + 2fi H a XiIi 



(3) 



where the first term is the energy of the circuit owing 
to the presence of the current distribution in the whole 
superconducting region, the second term is the energy 
due to the uniform applied magnetic fi eld , Mn are the 



self and mutual inductances from Eqs. (B5)-(B7), and 2, 



and Ij are the total current intensity that flows through 
the circuits labelled as i and j, respectively. Since no 
internal field dependence is considered for the critical 
current, \Ii\ = (Ax)(Az)J c . The sign of Ii is taken as 
positive when the current of the element at x > of the 
pair follows the positive y axis direction and negative 
otherwise. 

In the initial magnetization curve, after using the en- 
ergy minimization procedure for a given applied field H a 
to find the current profile, we can calculate the magne- 
tization, the total magnetic field, and the magnetic field 
lines directly form the current distribution. 
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The magnetization, denned as the magnetic moment 
per unit volume, has only one non-zero component, M z , 
which can be calculated as 



M, 



in. 
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(4) 



where m z is the total magnetic moment of the set of 
strips. 

The two nonzero components of the total magnetic flux 
density B z and B x are calculated as the addition of all 
the closed circuit contributions, which can be calculated 
analytically integrating the Biot-Savart law |2J|. 

The magnetic flux lines are calculated as in |l(J us- 
ing that for translational symmetry the level curves of 
the y component of the vector potential, for the gauge 
V • A = 0, can be taken as the magnetic flux lines. The y 
component of the vector potential results from the con- 
tribution of all the elements pairs A y ^, whose analytical 
expressions are the Eqs. (B2)-(B4). 

We will consider the application of an external AC field 
H a = H m cos(wi). In this case we can calculate the imag- 
inary part of the AC susceptibility, x"i defined as 
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dff a M rcv (tf a ). 



(5) 
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The energy loss W per AC cycle of amplitude H a is re- 
lated to x" by the expression 26 ^] 



W = wHix" 



(G) 



III. CURRENT PENETRATION AND FIELD 
LINES 

The most important issue to study in the system of 
superconducting strips is the influence of magnetic cou- 
pling. It is well known that the superconductors tend to 
shield the magnetic field change not only in their interior 
but also in the space between two superconducting re- 
gions (as illustrated in the classical example of a hollow 
superconducting tube, for which the whole volume inside 
the tube and not only the tube walls are shielded). In 
our array of superconducting strips one should observe a 
related effect, because the current induced in the strips 
will tend to shield the effect of the magnetic field in their 
interior as well in the space between them. To study this 
effect we present in Fig. 2 calculations corresponding to a 
vertical array of three strips with b/a = 0.1 and different 
separations (h = 2a, h = 0.2a, and h = 0.02a, respec- 
tively) . The different profiles correspond to applied fields 
-ff a /-ff P cn =0.2, 0.4, 0.6, 0.8, and 1, where H peD is the field 
at which the array is fully penetrated by current (the cal- 
culation of such a field is discussed in Appendix 1). The 



results show unambiguously the strong influence of mag- 
netic coupling in the case that the separation is small, as 
illustrated in the case for h — 0.02a for which the cur- 
rent profiles are almost the same as if there were no gaps 
between the strips. On the other hand, the case h = 2a 
is already an example of a very weak magnetic coupling 
that result in a current penetration for each strip almost 
as if the two others were not present. 

To further study this point we present in Fig. 3 the 
field lines calculated for the three arrays of Fig. 2, for 
applied fields H a /H pen = 0.4. The left images show the 
total field and the right ones only the field created by 
the currents in the superconductor. It is clear that for 
h = 0.02a and even h = 0.2a, the applied field in the 
space between superconductors is basically shielded by 
them, in contrast with the case of h = 2a, for which 
the magnetic field is modified near each strip but not 
enough to make a significant contribution to the other 
two strips. Another way of seeing this effect is illustrated 
in the calculations of the field created by currents (right 
images in each figure) . One can observe that in all cases 
currents create a basically constant field in the spaces 
between strips. However, when the separation is large 
the total field lines created by the current are wrapped 
around each tape, whereas when the separation decreases 
up to h/a = 0.02 the field lines are hardly distinguishable 
from the case of the three tapes forming a single thicker 
one. 



IV. MAGNETIZATION 

In this section we analyze the magnetization of the 
arrays, calculated from the currents following Eq. (Q). 
The reverse and returning curve can be obtained from 
the initial one using Eqs. ([!]) and (|J). 

There are several important properties of the tapes 
arrays that can be understood from the magnetization 
results. In fig. 4 we plot the calculated magnetization 
M as function of the applied field H a for a set of three 
arrays, each with semisides ratio b/a of 0.01 and with 
different separations h/a = 0.02,0.2, and 2, respectively. 
We also plot the magnetization for a single tape with 
b/a = 0.01, another one with b/a = 0.03 correspond- 
ing to the case that the three tapes are one on top of 
each other. One can observe that the general trend is 
that the saturation magnetization remains the same for 
all cases, whereas the initial slope of the magnetization 
curve changes. The slope is largest in absolute value for 
the case of a single thin tape (b/a = 0.01), has intermedi- 
ate values for the three arrays (the largest slope is for the 
array with large separation and smallest for the one with 
the smallest one), and finally the case with b/a = 0.03, 
which, interestingly, is hardly distinguishable from the 
array with the smallest separation. The reason for such 
an enhancement of the initial slope is the demagnetizing 
effect associated with the large sample aspect ratio, ac- 
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cording to which the thinner the sample the larger the 
initial slope |l3|]lC| l. In the arrays for which the separa- 
tion between the tapes is small, the magnetic coupling 
increase, in agreement with the discussion in section III, 
so that the sample is behaving as having a larger thick- 
ness, and therefore, less demagnetizing effects and less 
initial slope. In order to study in more detail this effect, 
we have included in Fig. 4 the calculated magnetization 
of a single tape with b/a = 0.05; this tape corresponds to 
the array with h/a = 0.02 but as if the gaps between the 
tapes were filled by superconducting material as well. We 
can see that this case has the smallest slope, and there 
are large differences with respect to the case of the array 
of three tapes with separation h/a = 0.02. From these 
results we can conclude that, provided that the strips are 
close enough, the behavior of a z array is similar to that 
of a strip of thickness the sum of the superconducting re- 
gions of the z-array, and not the sum of the whole z-array 
including the gaps. 

Another interesting feature to study is the effect of 
the addition of more strips to the array. We compare in 
Fig. 5 the calculated M(H a ) curves for arrays with a fix 
distance (h = 0.2) and different number of tapes. We 
include in the figure the two analytically known limits of 
one infinitely thin strip ji| and for an infinite set of strips 
iP^ij . Results show a practical coincidence between the 
calculated results for a single tape with finite although 
small thickness and the results from the analytical for- 
mula for very thin strips [Q. With adding more tapes, 
the initial slope of the magnetization (defined as mag- 
netic moment divided by volume, so independent of the 
superconductor volume) gets smaller in absolute value. 
We find that even the case of 25 strips is significatively 
different from the Mawatari case for an infinite stack, so 
we can conclude that Mawatari's formula should be valid 
only for a very large number of tapes. In Fig. 6 we show 
the magnetization curves M(H a ) for arrays of tapes sep- 
arated a smaller distance (h/a = 2). It can be observed 
that now the differences in the slope are smaller than in 
the previous case, because there is less magnetic coupling 
among the tapes. 

V. AC LOSSES 

In this section we study the imaginary part of the 
AC susceptibility, \", calculated form the magnetization 
loops from Eq. (Q). \" 1S directly related to the power 
losses by Eq. (||). 

We present in Fig. 7 the calculated results for x" 
as function of the maximum applied AC field H a cor- 
responding to the magnetization curves of Fig. 4. All 
curves show a peak at some value of the applied field 
amplitude. It can be seen that the peak corresponding 
to the maximum in x" (and therefore a change in slope 
of the AC losses, since they are proportional to x" times 
H a as seen in Eq. (^J)) is shifting to higher fields and 



decreasing in magnitude with decreasing separation dis- 
tance. The reason for that can be obtained from the 
analogous shifting in the initial slope of the magnetiza- 
tion shown in Fig. 4. Since the AC losses are related to 
the area of the M(H a ) curve [Eqs. (JsJ) and (J|)] and the 
magnetization saturation is the same for all cases, the key 
factor for the losses behavior is the initial slope governed 
by the demagnetizing effects as discussed in Section IV. 

Results in Fig. 7 show that the x"(H & curve goes as 
-ff™ for low fields, with n ranging from 1.5 for a strip 
with b/a = 0.01, corresponding to the high h/a limit for 
an x-array, to 1.3 for a strip with b/a — 0.03, being the 
low h/a limit for arrays. These values lie between the 
known limiting values for infinite slabs (n — 1) |l[] and 
thin strips (n = 2) p9| . The corresponding power losses 
(related to x" by Eq. (|^); see inset) go therefore as -ff™, 
with n ranging from 3.5 to 3.3. For high fields, x" goes 
as -ff" 1 in all cases. 

In Figs. 8 and 9 we calculate the results for x" as func- 
tion of the maximum applied AC field H a with the goal 
of studying the effect of adding more tapes in the array. 
We use the same cases as in Figs. 5 and 6, correspond- 
ing to two different separations in the arrays (h/a = 2, 
and 0.2, respectively). Again, there is a close relation 
between the increase in the initial slope of the magneti- 
zation curves and the shifting of the position of the peak 
in x" t° higher fields. It is interesting to comment, how- 
ever, that the formulas provided by Mawatari for the case 
of an infinite array are not adequate to describe quanti- 
tatively the AC losses of a finite array even if the array 
consists of up to 25 tapes. 

It is clear from the calculations presented in this sec- 
tion that in order to reduce the magnetic contribution 
to the AC losses in a real tape made of superconduct- 
ing filaments with a given thickness one should increase 
the coupling of the filaments by decreasing the distance 
between them (see Fig. 7). Also, the results in Fig. 8 
and 9 show that the AC losses decrease with increasing 
the number of filaments while keeping the distance be- 
tween them, particularly in the case of a small separation 
(h/a = 0.2, Fig. 9). These results can be understood 
again as arising from the effect of demagnetizing fields. 
When the superconducting tapes are separated a small 
distance, the whole tape is acting basically as a single 
superconducting tape with a thickness equal to the sum 
of the thickness of the tapes, so the demagnetizing effects 
are less, the initial magnetization has a smaller slope (in 
absolute value) as discussed in Section V, and the area 
of the hysteresis loops is less for a given maximum ap- 
plied AC field, so that the AC losses are reduced. On 
the other hand, when the separation of the tapes gets 
larger, the magnetic coupling between them is less so 
that they act more like magnetically decoupled filaments 
with high aspect ratio and the demagnetizing effects act 
more, increasing the AC losses. It is interesting to notice 
that when we talk about the separation h being small 
one should understand it as small compared with the 
horizontal dimension 2a and not to the strip thickness 
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2b. Calculations of the same values of h/a for the case 
b/a = 0.1 (instead of b/a — 0.01 as in the results pre- 
sented above), for example, yield qualitatively the same 
effects for the same h/a values. 



VI. CONCLUSIONS 

We have presented a numerical model for calculating 
current penetration and field profiles, and magnetization 
and AC losses of an array of superconducting tapes. In 
this work we have analyzed the case of an array of ver- 
tically arranged superconducting strips. We have found 
that the demagnetizing effects have strong influences in 
the magnetic response of the tapes and the AC losses ap- 
pearing when an AC field is applied. We find that AC 
losses are reduced when decreasing the vertical separa- 
tion between filaments. When the vertical separation is 
small as compared with the filaments width, then the 
tape is behaving as a single filament with thickness the 
sum of the superconducting material. These results could 
be used as guides for designing actual superconducting 
tapes. Then, in order to optimize the losses, for fila- 
ments with a fixed aspect ratio, it is preferable to have a 
large number of them separated small distances so that 
there is a good magnetic coupling between them, as has 
been already experimentally found |2§| ]. The model can 
be applied to horizontal and matrix arrangements as is 
presented in the next paper. 
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APPENDIX A: FIELD OF FULL PENETRATION 

In this section we present a simple way to analytically 
calculate the field of full penetration of the array, H pcn , 
defined as the minimum applied field in the initial mag- 
netization curve for which current fills the whole of the 
superconducting region. The penetration field can be 
calculated in general as minus the field generated by the 
current distribution Hj in the last induced current point, 
where Hj = Hjz |l(],|27j]. So, both the current distribu- 
tion at the penetration field and the last induced current 
point position r m must be known to calculate H pen . 

In the geometry of a set of rectangular cross-section 
strips ordered as a z-stack we can distinguish two differ- 
ent cases, depending on whether the number of strips rif 
is odd or even. Although in both cases the current dis- 
tribution at iJpen is evident, it is not so for r m . The last 



induced current point can be found as where the field gen- 
erated by the currents is maximum in magnitude, since 
r m is the last point where external field is shielded. For 
the odd number of strips simply the position 

of the center of the central strip, although it is not so 
easy to determine when the number of strips is even. 

For a z-stack of an odd number of strips nf with di- 
mensions 2a and 2b in the x and z directions separated 
a distance h the penetration field is therefore 



H pen (a,b,h,n f ) 



2tt 



Fi (0,o,6) + 



2 J2 Fi{(?b + h)i,a,b) 



where F\(u,t, d) is defined as 



(n f odd), (Al) 



I u + d u — d 
l<i[u,t,d) — 2t< arctan arctan 



(u — d) In 
[u + d) In 



(u - df 



t 2 + {u- d) 2 

(u + d) 2 
t 2 + (u+ d) 2 



(A2) 



Eq. (Al) has been derived using the expression for the 
magnetic field created by a completely penetrated strip 
with uniform J c calculated by direct integration of the 
Biot-Savart law. The case nj = 1 reproduces the known 
result for the penetration field for a strip []l0|| . 

As mentioned above, when a z-stack have an even num- 
ber of strips we must find the last point where current is 
induced, at which the self field Hj is maximum in mag- 
nitude. Owing to the symmetry of the current fronts in 
the yz plane, this point will be on the z axis. Thus, only 
maximization of Hj z along the z axis is needed. Since 
the minimization of the field Hj z is different for every 
specific value of «/, we only report the result for nj — 2, 
which is the most important case, concerning the mag- 
netic coupling, for an even nf. Then, the last induced 
current point will be at a position 



z 2 = 



1 



-a 2 -h/3 



y/(a? + h(3) 2 + hfi{2a 2 + h 2 /2 + 2p 2 + hp) , (A3) 

being z m the z component of r m and [3 defined as (3 — 
2b + h/2. The penetration field for two strips is 

H pcn (a, b,nf = 2) = \F 1 (z m - b - h/2, a, b) + 
Ztt 

F l {z m +b + h/2,a,b)], (A4) 
where the function Fi(u,t,d) is defined in Eq. (|A2[). 
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APPENDIX B: CALCULATION OF 
INDUCTANCES 

In this appendix we calculate the self and mutual in- 
ductances used in Eq. (||). These inductances are calcu- 
lated for closed circuits of the shape of a pair of identical 
rectangular infinite prisms of dimensions 2a' x 26' car- 
rying uniform current density. The prisms are set sym- 
metrically to the zy plane, taking the y axis parallel to 
the infinite direction. The current of the prism set in the 
x > region is taken positive, while it is taken negative 
for the other. 

The self and mutual inductances are calculated from 
the magnetic energy using the equation Eq] 



M.J.I, = Wij 



dx I dzA y j(x, z)Ji — 

xt—a' J Zi—b' 

Xi+a' /'Zi+b' 

dx I dzA y j(x,z)Ji, (Bl) 
i—a' J Zi—b' 



where My is the mutual inductance per unit length of 
two closed circuits labelled as i and j respectively, Ii and 
Ij are the current intensity flowing through the circuits, 
Wij is the magnetic energy of the circuits, (xi,Zi) is the 
central position of the prism in the x > region of the 
i circuit, a' and b' the dimensions of the prisms in the x 
and z directions respectively, and A y j is the y component 
of the vector potential created by the circuit j taking the 
gauge V • A = 0. 

The vector potential A y j can be calculated by direct 
integration leading to 



Ay^j (x, z) 



2vr 



[F(x-x j> z)-F(x + x j ,z)}, (B2) 



where the function F(u,v) is defined as 

F(u, v) = f(a' - u, b'-v) + f(a' - u, b' + v) + 

f(a' + u, b'-v) + f(a' + u, b' + v), (B3) 



defining f(t, d) as 



f(t,d) = - 



td\n (t 2 +d 2 ) - 3td + 



d 



t arctan — Yt arctan — 



t 



t 



(B4) 



Taking into account that the current density in the 
pris ms is unifor m, can be deduced integrating Eq. 
® using Eqs. @-(p3), which yield 



M - - ^ 

13 16^a' 2 6' 2 



[G(a 



•^ii Uj Vi) 

-G(-Xj - Xi,yj -yi)], (B5) 



where the function G(u, v) is defined as 



G(u,v)= J2 ("!)* 



fc,Z,n,m— 1 

g(R(k,n)a' + u,R{l,m)b' + v), (B6) 

defining R(i,j) = (— l) 1 — (— l) 3 ', and the function g(t, d) 
as 

, 24 o , 2 t 4 dt 3 d td 3 t 

q{t, d) = — t"d" arctan arctan — + 

yy ' ' 48 96 6 t 6 d 



1 
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FIG. 1. Sketch of the array of superconducting strips. 

FIG. 2. Current profiles for an array of three superconduct- 
ing tapes of width 2a and height 2b separated a distance (a) 
h/a=2, (b) h/a=0.2, and (c) h/a=0.02. The profiles corre- 
spond to applied fields H a /H pen =0.2, 0.4, 0.6, 0.8, and 1, 
where H pEn is the penetration field of the array. For the sake 
of clarity, the separation, thickness and width of strips are 
not on scale. 

FIG. 3. Field lines corresponding to an applied field 
Ha./ .Hpen=0.4, where H pcn is the penetration field of the 
z-stack, for the stacks of Fig. 2. Right and left figures corre- 
spond to the total and self-field magnetic field lines, respec- 
tively. The distances are d/a — 2 (a,b), 0.2 (c,d), and 0.02 
(e,f). 

FIG. 4. Initial magnetization as a function of the applied 
field for, from left to right: a single tape with b/a = 0.01, 
an array of three tapes with b/a = 0.01 separated a dis- 
tance h/a = 2, the same array with a separation distance 
of h/a — 0.2, the same array with a separation distance of 
h/a = 0.02, a single tape with b/a — 0.03, and a single tape 
with b/a = 0.05. In the inset complete magnetization loops 
are plotted for the first and last cases. 

FIG. 5. Initial magnetization as a function of the applied 
field for, from left to right: a single tape with b/a = 0.01, 
the theoretical expression for thin strips (almost overlapped), 
a set of 3 tapes with b/a — 0.01 each, a set of 5 tapes with 
b/a = 0.01 each, a set of 9 tapes with b/a — 0.01 each, a set 
of 25 tapes with b/a = 0.01 each, and the expression given by 
Mawatari (Ref. 15) for an array of infinite number of tapes of 
dimensions b/a — 0.01. The separation distance between the 
tapes in all cases is h/a=0.2. 



FIG. 7. Imaginary part of the AC susceptibility x" as a 
function of the amplitude of the AC field _ff a corresponding 
to the curves of Fig. 4, in the same order (the case of a 
single tape with &/a=0.05 is not plotted in this figure). The 
corresponding power losses are shown in the inset. 

FIG. 8. Imaginary part of the AC susceptibility \" as a 
function of the amplitude of the AC field i7 a corresponding 
to the curves of Fig. 5, in the same order. 

FIG. 9. Imaginary part of the AC susceptibility x" as a 
function of the amplitude of the AC field i7 a corresponding 
to the curves of Fig. 6, in the same order. 



FIG. 6. Same as Fig. 5, but for a separation h/a=2. 



7 




a 




1 



1 



1 



1 



1 



-L 



0.00 0.01 0.02 0.03 



0.04 0.05 0.06 

H a /(J„a) 



0.07 0.08 0.09 0.10 




I I I I I I I I I I I I I I I I I I I I I 

0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.10 



H a /(aJ c ) 




0.000 0.005 0.010 0.015 0.020 0.025 



H a /(aJ c ) 




H ac /(aJc) 




H a c/(H) 




H a c/(aJ c ) 



